THE ASTROPHYSICAL JOURNAL LETTERS, 970:L29 (16pp), 2024 August 1 
€ 2024. The Author(s). Published by the American Astronomical Society. 


OPEN ACCESS 


https: / /doi.org/10.3847/2041-8213/ad566a 


CrossMark 


JWST Reveals CO Ice, Concentrated CO, Deposits, and Evidence for Carbonates 
Potentially Sourced from Ariel’s Interior 


Richard J. Cartwright ©, Bryan J. Holler*®, William M. Grundy**®, Stephen C. Tegler*®, Marc Neveu^?6», Ujjwal Raut’ ©, 
Christopher R. Glein’®, Tom A. Nordheim! ©, Joshua P. Emery*®, Julie C. Castillo-Rogez®®, Eric Quirico’ ®, 
Silvia Protopapa'°®, Chloe B. Beddingfield! ©, Matthew M. Hedman’! ©, Katherine de Kleer'*®, Riley A. DeColibus®®, 


Anastasia N. Morgan" 6, Ryan Wochner*®, Kevin P. Hand®®, Geronimo L. Villanueva? ®©, Sara Faggi ° ©, 


Noemi Pinilla-Alonso'*®, David E. Trilling^ 9, and Michael M. Mueller? © 
! Johns Hopkins University Applied Physics Laboratory, 11101 Johns Hopkins Rd., Laurel, MD 20723, USA; richard.cartwright @jhuapl.edu 
T Space, Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA 
~ Lowell Observatory, 1400 W. Mars Hill Road, Flagstaff, AZ 86001, USA 
4 Department of Astronomy and Planetary Science, Northern Arizona University, 527 S. Beaver Street, Flagstaff, AZ 86011, USA 
~ University of Maryland, 4296 Stadium Dr., College Park, MD 20742, USA 
5 Solar System Exploration Division, NASA Goddard Space Flight Center, 8800 Greenbelt Road, Greenbelt, MD 20771, USA 
Z Southwest Research Institute, 6220 Culebra Road, San Antonio, TX 78238-5166, USA 
8 Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Drive, Pasadena, CA 91109, USA 
a Institut de Planétologie et d’Astrophysique (IPAG), UMR 5274, CNRS, Université Grenoble Alpes, Grenoble, France 
1 Southwest Research Institute, 1301 Walnut Street, Boulder, CO 80302, USA 
H Department of Physics, University of Idaho, 375 Perimeter Drive, MS 0903, Moscow, ID 83843, USA 
12 Division of Geological and Planetary Sciences, California Institute of Technology, 1200 E. California Boulevard, Pasadena, CA 91125, USA 
13 American University, 4400 Massachusetts Avenue NW, Washington, DC 20016, USA 
14 Florida Space Institute, University of Central Florida, 12354 Research Pkwy., Orlando, FL 32826, USA 
5m Kapteyn Astronomical Institute Rijksuniversiteit Groningen, P.O. Box 800, 9700 AV Groningen, The Netherlands 
Received 2024 May 7; revised 2024 June 8; accepted 2024 June 9; published 2024 July 24 


Abstract 


The Uranian moon Ariel exhibits a diversity of geologically young landforms, with a surface composition rich in 
CO, ice. The origin of CO, and other species, however, remains uncertain. We report observations of Ariel's 
leading and trailing hemispheres, collected with NIRSpec (2.87—5.10 um) on the James Webb Space Telescope. 
These data shed new light on Ariel's spectral properties, revealing a double-lobed CO; ice scattering peak centered 
near 4.20 and 4.25 jum, with the 4.25 um lobe possibly representing the largest CO, Fresnel peak yet observed in 
the solar system. A prominent 4.38 jum '*CO, ice feature is also present, as is a 4.90 um band that results from 
'2CO, ice. The spectra reveal a 4.67 jum "CO ice band and a broad 4.02 jum band that might result from carbonate 
minerals. The data confirm that features associated with CO, and CO are notably stronger on Ariel’s trailing 
hemisphere compared to its leading hemisphere. We compared the detected CO, features to synthetic spectra of 
CO, ice and mixtures of CO, with CO, H20, and amorphous carbon, finding that CO» could be concentrated in 
deposits thicker than ~10 mm on Ariel's trailing hemisphere. Comparison to laboratory data indicates that CO is 
likely mixed with CO». The evidence for thick CO; ice deposits and the possible presence of carbonates on both 
hemispheres suggests that some carbon oxides could be sourced from Ariel's interior, with their surface 
distributions modified by charged particle bombardment, sublimation, and seasonal migration of CO and CO; from 
high to low latitudes. 


Unified Astronomy Thesaurus concepts: Carbon dioxide (196); Ice spectroscopy (2250); Surface ices (2117); 


Surface processes (2116); Surface composition (2115); James Webb Space Telescope (2291); Uranian 


satellites (1750) 


1. Introduction 


Ariel is a candidate ocean world (e.g., Hendrix et al. 2019; 
Castillo-Rogez et al. 2023), with a young surface («1 Ga in 
some regions; Kirchoff et al. 2022), exhibiting large chasmata 
and putative cryovolcanic features (e.g., Smith et al. 1986; 
Beddingfield & Cartwright 2021; Beddingfield et al. 2022). 
Ground-based telescope studies determined that the surfaces of 
Ariel and the other large Uranian moons are composed of H5O 
ice mixed with a neutral absorber analogous to amorphous 
carbon (e.g., Cruikshank et al. 1977; Brown & 
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Cruikshank 1983; Clark & Lucey 1984), overprinted by 
deposits of "pure" CO, ice (Grundy et al. 2003, 2006; 
Cartwright et al. 2015, 2022), and possible ammonia (NH3) 
bearing compounds (e.g., Bauer et al. 2002; Cartwright et al. 
2018, 2020a, 2023; DeColibus et al. 2022, 2023). 

Because of the high obliquity of the Uranus system (~98°), 
the subsolar point migrates between 82? south and 82? north 
over the course of Uranus's orbit, exposing the summer poles 
of its moons to constant sunlight (estimated peak temperatures 
80—90 K; Hanel et al. 1986; Sori et al. 2017) while blocking 
sunlight from reaching their winter poles (20—30 K; Sori et al. 
2017), for about 21 yr each season. At the estimated peak 
temperatures during spring and summer, CO; can sublimate 
and likely migrates to the winter hemisphere. Integrating this 
process over many seasons, the predicted outcome is to strip 
the poles of CO, and other volatiles and concentrate them at 
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low latitudes, where diurnal variations in heating reduce 
sublimation rates (Grundy et al. 2006; Sori et al. 2017; 
Steckloff et al. 2022; Menten et al. 2023). Over time, CO, 
should be gradually lost to space owing to Jeans escape and 
magnetospheric interactions, depleting surface deposits, poten- 
tially enriching heavy isotopes, and raising the possibility that 
CO; is actively replenished. 

By analogy to surface irradiation of the Galilean satellites 
and Saturn's midsized moons, radiolytic generation of CO, and 
other species might occur on the Uranian satellites (Grundy 
et al. 2003, 2006; Cartwright et al. 2015). Such a process might 
explain the stronger CO, "triplet" band (1.9—2.1 um) on the 
trailing hemispheres of the largest Uranian moons and its 
weakening with increasing orbital distance (i.e., strongest on 
Ariel, weakest on Oberon). Nevertheless, predictions of moon- 
magnetosphere interactions at Uranus are not well constrained, 
and perhaps CO, is native and sourced from these moons’ 
interiors, with larger deposits on their trailing sides due to 
enhancement by radiolysis, similar to the multiple origin 
scenarios proposed for CO, on Jupiter's moon Callisto (e.g., 
Hibbitts et al. 2000, 2002; J. M. Moore et al. 2004; Cartwright 
et al. 2024). 

Irradiation of CO, should also generate other species, 
including CO (e.g., Bennett et al. 2010b; Raut et al. 2012; 
Raut & Baragiola 2013; Mifsud et al. 2022), carbon suboxide 
(C303; e.g., Strazzulla et al. 2007), and perhaps the cyanate ion 
(OCN ) and other CN-bearing compounds, assuming that CO; 
can interact with NH3 and H5O (e.g., Hudson et al. 2001). 
Supporting this scenario, some ground-based spectra of Ariel 
exhibit a 2.35 jum band, tentatively attributed to CO ice 
(Cartwright et al. 2022). While hitherto undetected, hydrogen 
peroxide (H5O;) may be present on Ariel, forming from 
irradiation of H2O, similar to Europa (e.g., Carlson et al. 1999) 
and Charon (Protopapa et al. 2023, 2024). Laboratory 
experiments demonstrate that H2O, formation is more efficient 
in H50 ice mixed with a small fraction of CO; (<9%; Mamo 
et al. 2023), potentially making the Uranian moons' surfaces 
ideal for H5O» production. 

Unlike the CO, triplet band, subtle features between 2.12 
and 2.27 um attributed to NH3-bearing species and ammonium 
salts (hereafter referred to as the “2.2 wm band") do not display 
discernible hemispherical or orbital trends. Instead, their 
variable band strengths suggest association with local-scale 
geologic features and terrains that cannot be resolved in the 
available disk-integrated data sets (Cartwright et al. 
2020a, 2023; DeColibus et al. 2023). If NH3-bearing species 
are present, they are likely replenished over short timescales, 
due to the predicted rapid decomposition of NH3 by magneto- 
spheric charged particles (~10° yr at Miranda; Moore et al. 
2007). The short life span of NH; and the possible association 
between geologic features and the 2.2 um band on Ariel could 
result from recent exposure of NHs-rich deposits (Cartwright 
et al. 2020a). However, the ancient surface of the neighboring 
moon Umbriel (74.5 Ga; Kirchoff et al. 2022) displays 2.2 uum 
bands as well, raising the possibility that more refractory 
components are contributing, such as carbonates, phyllosili- 
cates, and organic residues (Cartwright et al. 2023). 

To investigate the nature and origin of CO; ice and 
determine whether NH3-bearing species and other constituents 
are present, we measured spectra of Ariel with the NIRSpec 
spectrograph on the James Webb Space Telescope (JWST). 
These observations span an important wavelength range for 


Cartwright et al. 


detecting CO; and CO (4.2-5.0 um), NH-bearing species 
(2.9-3.2 um), hydrocarbons (3.2-3.7 um), nitriles 
(4.3—4.8 um), carbonates (3.9—4.1 um), and H20, (~3.51 um). 


2. Data and Methods 


NIRSpec Observations: As part of General Observer (GO) 
Program 1786, NIRSpec on JWST (Jakobsen et al. 2022; 
Böker et al. 2023) observed Ariel with the G395M/F290LP 
grating (2.87-5.10 um, resolving power R ~ 1000) on 2023 
September 6 and 7, when the mid-observation, subobserver 
longitude was near 293° west (trailing hemisphere) and 63° 
west (leading hemisphere), respectively (subobserver latitude 
~64.°5 north for both observations). Each observation 
consisted of four dithers with NIRSpec’s integral field unit 
(IFU), for a total of ~3735 s (leading hemisphere) and ~3793 s 
(trailing hemisphere) of on-target time, using the NRSIRS2- 
RAPID readout mode. The data were downloaded from the 
Mikulski Archive for Space Telescopes (10.17909 /cwsn-7z50) 
(data reduction procedures described in Appendix A.1). 

Band Measurements: We measured the band areas and 
depths of detected absorption features and scattering peaks 
using a band measurement program that defines and divides off 
a local linear continuum for each feature (e.g., Cartwright et al. 
2024). Depth and height measurements were made by 
averaging the reflectance values within +0.002 to 0.003 jum 
of a user-defined band (B) or peak (P.) center, and 
uncertainties were computed using standard error propagation 
procedures (e.g., Taylor 1997). The spectral contrast for each 
absorption band (1 — B4) and scattering peak (P4 — 1) was then 
measured. The trapezoidal rule was used to measure the area of 
each feature, using Monte Carlo simulations sampling the 1c 
errors for data points within each feature to estimate errors 
(Table 1). 

Comparison to Synthetic and Laboratory Spectra: We 
compared Ariel's spectral properties to one-layer Hapke-Mie 
spectral models (Appendix A.2) generated using three sets of 
laboratory-derived indices of refraction (“optical constants") 
for crystalline CO» ice measured at 150 K (Hansen 1997), 70 K 
(Gerakines & Hudson 2020) and 21K (Quirico & 
Schmitt 1997a, 1997b) and for crystalline CO ice measured 
at 20K (Robert Brown, private communication; Gerakines 
et al. 2023). Particulate and areal mixtures of CO2, CO, H50, 
and amorphous C are shown in Figure A1. 

We present absorbance spectra of CO + CO» ice mixtures 
measured in the Astrophysical Materials Laboratory at North- 
ern Arizona University (Tegler et al. 2024). These co- 
condensed ices were made by mixing CO and CO, gases at 
room temperature before condensing as thin films on a gold 
mirror held at 20 K (see Tegler et al. 2024; Grundy et al. 2024a, 
for more details on laboratory procedures). We also compared 
the Ariel data to spectra of radiolyzed CO, ices processed by 
100 keV protons (Raut & Baragiola 2013). 


3. Results and Analyses 
3.1. Detected Spectral Features 


CO» Ice: The NIRSpec data exhibit a large scattering peak 
centered near 4.20 um on Ariel's leading and trailing hemi- 
spheres, flanked by an absorption band centered near 4.27 um, 
and the trailing hemisphere shows an additional peak near 
4.25 um (Figure 1). All of these features are associated with the 
asymmetric stretch fundamental (v3) mode of PROS (Figure 2). 
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Table 1 


Band Measurements for Ariel 


Cartwright et al. 


Feature Feature Wavelength >30 Spectral Contrast 
Feature Name Center Range Hemisphere Spectral Contrast Band Area and Band Area? Constituents 
(um) (um) (96) (107^ um) 
?3.10 jum 3.100 3.014-3.168 Leading 99.29 + 1.29 564.08 + 2.07 Yes H20 Ice 
Trailing 68.37 + 0.90 388.65 + 1.68 Yes 
3.33 um 3.331 3.308-3.362 Leading —0.65 + 0.68 —2.54 + 0.63 No "CO, Ice, CH4, 
Trailing 2.59 + 0.34 5.69 + 0.50 Yes Hydrocarbons 
4.02 um 4.021 3.954—4.109 Leading 3.80 + 0.31 35.47 + 1.17 Yes "CO, Ice, CO; 
Trailing 7.30 + 0.66 65.38 + 1.14 Yes 
74.15 um 14.147 4.117-4.187 Leading 7.33 + 0.48 35.90 + 0.86 Yes "CO, Ice, O-D 
74.150 4.117-4.178 Trailing 9.50 + 0.77 35.73 + 0.67 Yes 
*^4 20 uum 4.202 4.152-4.276 Leading 38.47 + 0.83 193.14 + 1.61 Yes CO, Ice 
Trailing 126.13 + 1.30 743.91 + 2.19 Yes 
425 um 4.252 4.152-4.276 Leading 15.64 + 0.85 193.14 + 1.61 Yes CO, Ice 
Trailing 82.42 + 1.10 743.91 + 2.19 Yes 
74.27 um 4.273 4.251-4.303 Leading 25.86 + 1.21 43.45 + 0.91 Yes CO, Ice 
Trailing 41.24 + 0.79 89.45 + 0.79 Yes 
“4.30 um 4.298 4.2924.302 Leading 3.22 + 1.43 1.35 + 0.52 No 60"C!80 Ice, 
Trailing 4.13 + 0.87 1.59 + 0.34 Yes amorphous CO2 
4.38 um 4.379 4.365-4.391 Leading 7.82 + 0.93 10.81 + 0.68 Yes CO, Ice 
Trailing 11.12 + 1.04 16.19 + 0.61 Yes 
441 um “4.403 4.395-4.423 Leading 2.96 + 0.88 3.20 + 0.56 Yes 15930180 Ice, 
74412 4.393-4.429 Trailing 5.26 + 1.13 8.38 + 0.71 Yes CxO», Nitriles 
4.47 um 4.465 4.430-4.495 Leading 2.66 + 0.86 10.22 + 1.05 Yes "CO, Ice, CxO», 
Trailing 6.11 + 0.63 27.06 + 1.04 Yes Nitriles 
4.59 um 4.593 4.5524.639 Leading 1.53 + 0.60 6.08 + 1.24 No "CO, Ice, OCN- 
Trailing 2.50 + 0.73 14.90 + 1.02 Yes 
4.67 um 4.675 4.640—4.706 Leading 6.60 + 1.09 11.64 + 1.03 Yes CO Ice 
Trailing 23.29 + 0.76 63.06 + 0.93 Yes 
4.78 um 4.780 4.764—4.795 Leading 1.70 + 0.65 0.35 + 0.57 No CQ Ice, HCN 
Trailing 2.42 + 0.25 4.17 + 0.78 Yes 
4.90 um 4.898 4.845-4.927 Leading 4.59 + 0.76 11.57 + 1.30 Yes 2CO, Ice, 
Trailing 20.29 + 0.96 64.77 + 1.06 Yes CO3 


Notes. L = Leading, T = Trailing. In the last column, confirmed constituents are shown in bold, and suggested constituents are italicized. 


* Spectral features are peaks measured above the continuum. 


P Band areas for the 4.20 and 4.25 jum features are convolved, and we report the same band area for both peaks. 

© The 4.15 um and 4.27 um features overlap the wavelength range of the 4.20 jum and 4.25 jum scattering peaks, and the reported measurements for these features (and 
their hemispherical asymmetries) are provided for completeness but are likely inaccurate. 

d The 4.30 jum feature is partly embedded in the long-wavelength wing of the 4.27 jum feature. 


Ariel’s leading and trailing sides also exhibit prominent 
absorption features near 4.90 um, which corresponds to a 
biphonon 4- phonon combination mode, resulting from 
collective vibrations across a CO, ice lattice (Bini et al. 
1991). Radiolytic carbon trioxide (CO3) may provide minor 
contributions to the 4.90 jum feature (Figure A2). Another 
absorption band centered near 4.38 um is consistent with the v4 
mode of "CO; ice. Ariel's large 4.25 um scattering peak 
coincides with the wavelength range where the extinction 
coefficient for CO», k, is >1 and the refractive index, n, is «1 
(Gerakines & Hudson 2020; Figure A3), suggesting that it is 
the largest CO» ice Fresnel peak yet observed on an icy body. 
For comparison, JWST revealed a much smaller CO, Fresnel 


peak on Charon, shifted to ~4.265 jum (Protopapa et al. 
2023, 2024). 

We calculated that 26.796 of Ariel's disk is sampled by both 
the leading and trailing hemisphere observations (Holler et al. 
2016), with most of this overlap (18.7%) at >45° north, where 
CO; is likely depleted (Grundy et al. 2006; Sori et al. 2017; 
Steckloff et al. 2022; Menten et al. 2023). We also measured a 
suite of less prominent absorption bands that may result from 
CO», C302, carbonates, and nitriles (Appendix A.7), centered 
near 3.33, 4.02, 4.15, 4.30, 4.41, 4.47, and 4.59 um (Figures 2 
and A4). Other tentative features that may result from CO» ice, 
but are too subtle to reliably measure, are centered near 3.01, 
4.84, and 4.93 um (Figures 2, 4, and A4). 
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Figure 1. Top: NIRSpec IFU (G395M/F290LP) reflectance spectra and 1c uncertainties for Ariel's leading (purple) and trailing (red) hemispheres, vertically offset 
for clarity and normalized to 1 at 4.145 jum. Bottom: close-up of the same spectra, focusing on the 3.8-5.1 jum wavelength region. Confirmed absorption bands and 
scattering peaks associated with H5O ice, CO, ice, and CO ice are shown in bold, whereas suggested constituents are italicized and followed by a question mark. 
Amorphous CO, is abbreviated to “A. CO;." The central wavelengths (jum) for the identified features are listed vertically along each dotted line and in Table 1. 


CO Ice: The NIRSpec spectra show a prominent absorption 
feature near 4.67 jum that is consistent with the v, mode of ^CO 
ice (e.g., Sandford et al. 1988; Bennett et al. 2010b), observed 
on a wide variety of trans-Neptunian objects (TNOs) and 
Centaurs by NIRSpec (e.g., Brown & Fraser 2023; Licandro 
et al. 2023; Emery et al. 2024; de Pra et al. 2024; Souza- 
Feliciano et al. 2024). CO ice should sublimate rapidly at 
Ariel’s peak surface temperatures (80-90 K), and it is probably 
replenished on short timescales and perhaps complexed with a 
less volatile component, such as CO» (e.g., Sandford et al. 
1988). Furthermore, Ariel's 4.67 jum feature exhibits flanking 
sidebands that are absent from pure CO ice (Figure 3, 
Section 3.3). A subtle feature near 4.78 jum on Ariel’s trailing 
hemisphere likely results from CO (see Appendix A.7 for 


other interpretations). We also report CO / "CO isotopic 
ratios (Appendix A.9) and modeling of possible "°C enrich- 
ment (Appendix A.10). 

H30 Ice: The spectral signature of H5O ice is apparent on 
both sides of Ariel, with clear evidence for the strong 3.0 um v4 
mode, the 3.1 jum Fresnel peak indicative of crystalline H20 
ice, the 3.6 wm H5O ice continuum peak, and a broad 4.5 um v; 
+ vg combination mode (Figure 1; Mastrapa et al. 2009 and 
references therein). 

What about NH3, Hydrocarbons, and H202? The NIRSpec 
spectra do not display evidence for the 2.96 um v3 mode 
indicative of NH-bearing species, nor other features associated 
with NH3 or NH4. NH-bearing features between 2.9 and 
3.2 um are also apparently absent from Charon (Protopapa 
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Figure 2. Comparison between Ariel’s leading (purple) and trailing (red) hemisphere spectra (vertically offset for clarity, normalized to unity at 4.145 um) and lo 
uncertainties, along with synthetic spectra of crystalline CO; ice (1-3) and CO ice (4). These synthetic spectra were generated using CO» grain diameters and optical 
constants from: model 1, 10 jum, Hansen (1997); model 2, 10 jum, Quirico & Schmitt (19972, 1997b); model 3, 50 jum, Gerakines & Hudson (2020); model 4, 10 jum, 
Robert Brown (private communication). Models 1 and 2 were scaled to the depth of the 4.90 jum band on Ariel’s trailing side, Model 3 was scaled to the height of the 
4.20 um scattering peak on Ariel’s trailing side, and Model 4 was scaled to the depth of the 4.67 jum band on Ariel’s trailing side. All four models are offset vertically 
for clarity. The central wavelength of measured features (Table 1) and possible features near 3.01, 4.84, and 4.93 um are indicated by dotted lines. 


et al. 2023, 2024), which exhibits a prominent 2.2 um band 
long attributed to NH-bearing compounds (e.g., Brown & 
Calvin 2000; Buie & Grundy 2000; Grundy et al. 2016; Cook 
et al. 2018, 2023; Protopapa et al. 2020). Thus, NH-bearing 
species may contribute to Ariel's 2.2 jum band but are obscured 
by H50 ice absorption in the 3 jum region sampled by G395M 
data. We also find no reliable evidence for C-H stretching 
modes exhibited by hydrocarbons between 3.2 and 3.5 um, nor 
any evidence for a prominent H2O, combination mode (v, + 
ve) near 3.51 um (Figure A5; Bain & Giguére 1955). 


3.2. Band Parameter Measurements 


We measured seven confirmed features in both spectra and 
another eight likely features on Ariel's trailing hemisphere, 
four of which are also on Ariel's leading side (>3o detection; 
Table 1). The 3.1 im H5O Fresnel peak is significantly stronger 
on Ariel's leading side. In contrast, all of the features that are 
definitively associated with CO, and CO ice are stronger on 
Ariel's trailing hemisphere. The distributions of H5O and CO; 
are consistent with the hemispherical trends established in 
ground-based studies (e.g., Grundy et al. 2003; Cartwright et al. 
2022). The features centered near 4.02, 4.41, and 4.47 um are 
significantly stronger on Ariel’s trailing side (3c difference). 
The 3.33, 4.30, 4.59, and 4.78 um bands are only reliably 
detected on Ariel’s trailing hemisphere (»3o detection; 
Table 1). 


3.3. Synthetic and Laboratory Spectra of CO» and CO Ice 


We compared the Ariel data to synthetic spectra of 
crystalline CO» ice (models 1-3) and CO ice (model 4), 
shown in Figure 2. These comparisons highlight the large 


number of subtle features that can be expressed by CO; ice 
(models 1 and 2), many of which appear to be present on Ariel 
(Table 1). Furthermore, Ariel's CO» scattering peaks can be 
approximated by a layer of large CO, grains (250 jum 
diameters; model 3). The synthetic spectrum of CO ice matches 
the central wavelength positions of Ariel's 4.67 and 4.78 um 
features (model 4), confirming the presence of this molecule 
and possibly confirming its heavy-carbon isotopologue. 

Figure 3(a) shows laboratory data of CO and CO; mixtures, 
whereas Figure 3(b) shows CO generated in proton-irradiated 
CO, ice films containing CO and O;, with trace amounts of 
CO; and O, (Raut & Baragiola 2013). The pure CO sample 
exhibits a narrow sideband that likely results from a long- 
itudinal optical (LO) phonon mode caused by collective 
oscillations within the ice's structure. Similar LO modes have 
been reported for CO» samples illuminated at oblique angles 
(Cooke et al. 2016). However, this 4.658 um sideband 
exhibited by Ariel's 4.67 jum feature is broader than pure CO 
and more similar to the sidebands exhibited by mixed CO + 
CO, samples, consistent with prior work (Sandford et al. 1988). 
We speculate that these sidebands result from CO dimers and 
trimers with vibrational frequencies distinct from mono- 
meric CO. 

The v4 mode in the proton-generated CO sample we 
analyzed is centered at shorter wavelengths (724.668 um) than 
Ariel's 4.67 um band (4.675 um) and does not match its 
sidebands (Figure 3). Similarly, another study that generated 
CO via electron irradiation of CO, (20-50 K) measured a band 
center of 4.671 um for its v3 mode (Mifsud et al. 2022), slightly 
offset from Ariel's 4.67 um band. Thus, radiolytic CO may not 
be the primary contributor to CO on Ariel. Alternatively, ice 
annealing processes and CO diffusion to lower energy sites in 
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Figure 3. (a) Comparison between a continuum-subtracted version of the 
4.67 um absorption band on Ariel’s trailing hemisphere (black) and laboratory 
absorbance spectra of pure CO ice and CO + CO» ice mixtures measured in the 
Astrophysical Material Laboratory at temperatures between 20 and 30 K, with 
the CO feature persisting at temperatures up to 80 K in the mixed CO + CO, 
samples (Tegler et al. 2024; Grundy et al. 2024a). Pure CO; ice is essentially 
featureless in this wavelength range and omitted for clarity. The CO;:CO 
mixture (9:1, green, centered near 4.674 jum) provides a good match to the 
central wavelength of Ariel’s 4.67 4m band, whereas mixtures with 
comparable amounts of CO, and CO (6:4, blue, centered near 4.672 jum) 
provide less ideal matches. Pure CO ice (red, 4.677 jum) is much narrower, 
offset to longer wavelengths, and exhibits a narrow sideband (4.667 jum), likely 
resulting from a longitudinal optical phonon mode, which does not match the 
broad short-wavelength sideband of Ariel’s 4.67 jum band. Neither the pure CO 
nor CO 4- CO, mixtures are able to match Ariel's long-wavelength sideband 
centered near 4.69 um. (b) Comparison between a continuum-subtracted 
version of the 4.67 jum band on Ariel's trailing hemisphere (black) and 
laboratory spectra of a sample dominated by CO ice generated via irradiation of 
a CO, thin film with 100 KeV H* protons at 25 K (light blue) and 50K 
(orange), shown in optical depth units (—In(R/Ro), where R is the reflectance of 
the film and gold substrate and Rọ is the reflectance of the bare gold substrate; 
Raut & Baragiola 2013). All spectra are normalized to 1 at their maxima. The 
radiolytic CO features do not provide good matches to the central position, 
width, or sidebands of Ariel's 4.67 jum band. 


the surrounding CO; lattice might obscure the signature of 
radiolytic CO over time, making interpretation more difficult. 

The other sideband (74.687 um) expressed by Ariel’s 
4.67 um feature is absent from pure CO ice, CO + CO; 
mixtures, and radiolytic CO, suggesting mixing with other 
species, including H5O (Sandford et al. 1988), CO clathrates 
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(Dartois et al. 2011), or perhaps carbonates. Comparison 
between these Ariel data and laboratory spectra of irradiated 
carbonates and other carbon oxides could test our inference, 
based on Figure 3, that CO may be (partially) native to Ariel. 
Similarly, new optical constants for CO + CO» mixtures and 
CO, mixed with other carbon oxides are likely required to 
better understand Ariel's spectral properties. 


4. Discussion 
4.1. Concentrated Deposits of CO; Ice 


The results presented here and in prior studies (Grundy et al. 
2003, 2006; Cartwright et al. 2015, 2020b) indicate that CO; 
on Ariel has remarkably similar spectral properties to crystal- 
line CO, ice measured in the laboratory (i.e., CO; molecules 
primarily bonded to each other in a long-range order). 
Comparison to synthetic spectra of crystalline CO, ice 
(Figure 2) highlights that Ariel exhibits many subtle CO; 
features outside of the strong v4 mode. Specifically, the 3.33 
jam and 4.90 um bands and other possible features near 3.01, 
3.33, 4.84, 4.90, and 4.93 jum on Ariel's trailing hemisphere do 
not appear on many other objects where CO; has been detected 
by JWST/NIRSpec (e.g., Brown & Fraser 2023; Pinto et al. 
2023; Villanueva et al. 2023; Bockelee-Morvan et al. 2024; 
Cartwright et al. 2024; de Pra et al. 2024; Emery et al. 2024; 
Protopapa et al. 2024; Wong et al. 2024; TNO examples 
provided in Figure A6), with the exception of Triton (Wong 
et al. 2023). Consequently, Ariel's surface exhibits some of the 
most CO»-rich deposits in the solar system. 

To estimate the thickness of Ariel's CO, ice deposits, we 
calculated the e-folding depths (1/a, where o is the absorption 
coefficient) for near-infrared (NIR) photons propagating 
through a slab of CO; ice, which is essentially an application 
of the Beer-Lambert absorption law. To account for scattering 
off grain boundaries, we calculated the mean optical path 
length, MOPL — —1/(a"In(R), where R is the reflectance 
(Clark & Roush 1984), using synthetic spectra composed of 
CO, ice grains (1, 10, and 100 jum diameters). The results of 
these two approaches indicate that photons spanning the 
wavelength range of the CO; v, band (4.1—4.4 jum) can travel 
0.0001—0.1 mm into slabs of CO, ice or deposits dominated by 
CO, grains before being absorbed (Figure 4). In contrast, 
photons penetrate 0.3 mm in the wavelength range of the 
4.90 um band, ~4 mm in the wavelength ranges of the 3.33 and 
4.84 um bands, and upward of 10 mm in the wavelength ranges 
of the possible 3.01 and 4.93 wm features (Figure 4). The 
penetration depth estimates reported here are broadly consistent 
with prior estimates made on weak CO; features between 
1.5-1.7 um (4—30 mm depths) and 1.9-2.1 um (0.07-0.4 mm 
depths) detected on Ariel (Cartwright et al. 2015, 2022). 
Additionally, these estimates are consistent with laboratory 
measurements of CO, ice optical constants, which used thin 
films (0.25—3 um thick) to avoid saturation while measuring the 
v4 mode (Gerakines & Hudson 2020) and samples >10 mm 
thick to measure weaker CO» features (Hansen 1997, 2005). 

As an additional test, we generated synthetic spectra 
composed of CO; in areal and particulate mixtures with CO, 
H2O, and amorphous carbon. We found that the CO, v3 mode 
is present in all areal and particulate mixtures, and weaker CO; 
features are present in areal mixtures (see Appendix A.2 for 
model details). In contrast, weak features near 3.01, 3.33, 4.84, 
and 4.93 um are not present in any of the particulate mixture 
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Figure 4. The e-folding photon penetration depths into a slab of CO; ice (red: Gerakines & Hudson 2020; purple: Quirico & Schmitt 1997a, 1997b), CO ice (orange; 
Gerakines et al. 2023), and H5O ice (black; Mastrapa et al. 2009), and MOPL for photons into a layer of CO» ice grains with 1 jum diameters (blue). The MOPLs for 
photons into layers composed of 10 and 100 jum diameter CO) ice grains are intermediate between the blue and purple lines. The black dashed lines highlight the band 
centers for Ariel's 3.33 and 4.90 jum bands, and the gray dashed lines highlight other subtle CO; ice features that may be present in the spectrum of Ariel’s trailing 
hemisphere (centered near 3.01, 4.84, and 4.93 um). This plot shows that only the strong v3 modes for ?CO,, PCO,, "CO, and '^CO should be exhibited by Ariel, if 
carbon oxides are well mixed with HO on its surface. The presence of weaker CO, ice features supports the presence of concentrated CO; ice deposits exposed on 
Ariel's surface (likely mixed with small amounts of CO and perhaps other carbon oxides). 


models, while the 4.90 um band is present in particulate 
mixture models that include >40% CO, (Figure Al). These 
results support the hypothesis that CO; ice on Ariel's surface is 
segregated from HO ice. Whether these CO,-dominated 
deposits are spatially associated with specific geologic land- 
forms, or are more regionally dispersed, cannot be determined 
with these disk-integrated spectra. 


4.2. Radiolytic Production and Seasonal Migration of CO and 
CO» 


At winter temperatures of 20—30 K, CO; and CO ices are 
stable and likely form a winter cap. Once exposed to sunlight in 
spring, CO, and CO should start to sublimate and migrate to 
cold traps near the new winter pole or leave the surface 
environment via Jeans escape, in particular for the more 
volatile CO. Seasonal migration of CO; molecules could lead 
to the formation of a transient layer of pure CO» frost grains 
that is predicted to be up to 2.4 mm thick (Steckloff et al. 
2022). Scattering within such a thick layer of seasonally mobile 
CO, grains might contribute to Ariel's large 4.20 and 4.25 um 
scattering peaks. Furthermore, the JWST/NIRSpec results 
presented here indicate that CO, ice deposits could be upward 
of 10 mm thick on Ariel's trailing hemisphere (Section 4.1), 
raising the possibility that CO» deposits on its trailing side are 
sufficiently thick to be partly resistant to seasonal migration. 
On Ariel’s leading hemisphere, CO, deposits may be only 
~0.3 mm thick (Figure 4), suggesting that CO» on its leading 
side may primarily result from sublimation and transport from 
its trailing hemisphere. 

One way to explain how Ariel retains CO is if it forms at 
depth via radiolytic decomposition of thick CO» deposits. 
Measurements made by Voyager 2 indicate that heavy ions are 


largely absent from Uranus's magnetosphere (Ness et al. 1986), 
suggesting that high-energy electrons and protons could be the 
primary drivers of radiolytic processes. Energy deposition by 
protons and heavy ions is mostly limited to the top 0.01 mm of 
icy regoliths (e.g., Delitsky & Lane 1998), but energetic 
electrons (~1 MeV) can penetrate centimeter-scale depths into 
ices (e.g., Nordheim et al. 2017). CO molecules generated at 
depth would be stable at Ariel's winter pole before diffusing 
out of its regolith once exposed to sunlight, possibly being 
retained long enough to be detected. Alternatively, geologic 
sources of native CO could help replenish CO ice on Ariel's 
surface (Figure 3, Section 3.3). 


4.3. Possible Internally Derived Materials 


The 2.2 jum bands detected in ground-based observations of 
the Uranian moons have been attributed to potentially 
internally derived deposits rich in NH; hydrates and NH; ice, 
carbonates, phyllosilicates, and/or organics (e.g., Cartwright 
et al. 2020a, 2023; DeColibus et al. 2023). Although we did not 
detect evidence of NH-bearing species, phyllosilicates, or 
hydrocarbons in the G395M data, Ariel’s 4.02 um feature 
could result from a v, 4- v4 combination mode expressed by 
CO; in carbonate minerals (Hexter 1958; Bishop et al. 2021), 
similar to Callisto's 4.02 m band (Johnson et al. 2004; 
Cartwright et al. 2024). Laboratory experiments show that 
although radiolytic CO; exhibits a 4.89 um feature (Figure A2; 
Raut & Baragiola 2013), a complementary 4 um CO; band was 
not observed in these irradiation experiments. Similarly, 
carbonic acid (H52CO3), generated via irradiation of H5O and 
CO; mixtures at cryogenic temperatures («100 K), exhibits a 
broad absorption band between 3.8 and 3.9 um (e.g., Moore & 
Khanna 1991; Hage et al. 1998; Gerakines et al. 2000), but it is 
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Figure A1. Synthetic spectra (Appendix A.2) composed of areal (models 1-3) and particulate (models 4—8) mixtures, offset vertically for clarity. These spectral 
models include crystalline H5O ice (Mastrapa et al. 2009), crystalline CO% ice (Quirico & Schmitt 1997a, 1997b), crystalline CO ice (Robert Brown, private 
communication), and amorphous carbon (Rouleau & Martin 1991), with grain diameters of 2, 10, 10, and 1 jum, respectively. Each model includes 1096 CO ice and 
0.596 amorphous carbon mixed with (1) 84.5% H5O ice and 5% CO; ice, (2) 79.5% H,O ice and 10% CO, ice, (3) 64.5% H,O ice and 25% CO; ice, (4) 64.5% H,O 
ice and 2596 CO» ice, (5) 49.596 H20 ice and 40% CO; ice, (6) 39.5% H20 ice and 5096 CO» ice, (7) 14.5% H20 ice and 75% CO» ice, and (8) 14.5% H5O ice and 
75% CO, ice. Model 7 includes 25% CO, (Quirico & Schmitt 1997a; 10 jum grains) and 50% CO; (Gerakines & Hudson 2020; 100 jum grains). These models 
demonstrate that the Ariel features between 4.19 and 4.30 jum and between 4.65 and 4.80 um, which are associated with the v3 modes of CO; and CO ice, 
respectively, are exhibited by all eight models (dotted lines with bold text). In contrast, weak CO» features tentatively observed on Ariel, near 3.01, 3.33, 4.84, and 
4.93 um (dotted lines, italicized text), are only exhibited by areal mixtures that include 21096 CO; (models 2-3) and are not observed in any of the particulate 
mixtures (2596—7596 CO» ice; models 4-8). The 4.90 jum band is expressed by all areal mixtures (596-2596 CO», models 1—3) and particulate mixtures with 24096 
CO, (models 5-8). The presence of weak CO, ice features in the areal mixtures and their absence from the particulate models demonstrate that concentrated deposits 
of pure CO, on Ariel’s surface are required for these features to be expressed (Section 4.1). None of the areal or particulate mixtures exhibit a ~4 um feature, 


suggesting that carbonates, or another component, are more likely to contribute to Ariel’s 4.02 jum band compared to CO; ice (dashed line, italicized text). 


uncertain what processes might cause this feature to shift to 
longer wavelengths, matching Ariel's 4.02 um band. Radi- 
olytically-formed, CO3-bearing species are therefore unlikely 
to contribute to Ariel’s 4.02 jum band, and nonradiolytic, native 
carbonates seem to be more likely to contribute. 

A prior study suggested that Umbriel's 2.2 um band may 
result from thermonatrite (Na?CO; : H20), possibly contribut- 
ing to bright crater floor deposits (Cartwright et al. 2023), along 
with cold-trapped CO; ice (Sori et al. 2017). In this scenario, 
emplaced carbonates could serve as base material for radiolytic 
CO2, or perhaps endogenic CO, could be delivered from 
Ariel's interior with carbonates, either at present or in the 
geologic past. Indeed, it is predicted that within deep oceans 
evolved from carbon-rich ices, a large fraction of CO, (several 
hundreds of mmol/(kg H5O)) may be in solution with 
bicarbonate (Castillo-Rogez et al. 2022). If CO; is outgassed 
from Ariel’s interior, then some of it should condense and 
contribute to the concentrated CO, ice deposits detected in 
JWST and ground-based data sets. Endogenic CO», if mixed 
with liquid H5O and coexisting with bicarbonate, could 
indicate a mildly acidic ocean (pH ~ 6-8), favoring the surface 
precipitation of bicarbonate salts and potentially Mg or Ca 
carbonates at lower pH values (e.g., Glein et al. 2015; Tosi 
et al. 2024). Alternatively, radiolytic CO», or native CO, 
sequestered in Ariel’s crust, would permit a more basic ocean 


chemistry (pH > 9-10), supporting Na and NH, carbonate 
precipitation. 

The presence of carbonates would have important implica- 
tions for minerals possibly formed in an aqueous environment 
and for the habitability of Ariel's interior, including the 
availability of phosphorus (Postberg et al. 2023), a key 
chemical component for life. Except for NH4 carbonates, 
carbonate salts require environments where silicate minerals 
can interact with liquid H5O (Castillo-Rogez et al. 2022). Such 
formation environments may exist at rock—water interfaces in 
ocean world interiors such as Ceres, whose strong 4 um band 
results from carbonates (e.g., Rivkin et al. 2006), including 
NaCO; salts (e.g., De Sanctis et al. 2016; Carrozzo et al. 2018; 
Raponi et al. 2019), likely formed from the alteration of rock in 
contact with an ocean that included NH3 (Castillo-Rogez et al. 
2018). Similarly, Enceladus's plume particles are dominated by 
H50 ice mixed with minor amounts of Na carbonates (Postberg 
et al. 2009) and Na phosphates (Postberg et al. 2023) that likely 
formed in its ocean. 

If CO is internally derived and released during winter, it 
could persist on Ariel’s surface, condensing with CO». To 
survive in Ariel’s interior, endogenic CO would need to be 
sequestered in its crust and unable to interact with liquid H20, 
or else it would have oxidized to CO» or reduced to metastable 
organic compounds, such as formic acid/formate and perhaps 
CH, (Neveu et al. 2015; Glein & Waite 2018). CO trapped as 
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Figure A2. Comparison between a continuum-subtracted version of Ariel's 
4.90 um band on its trailing hemisphere and a 4.89 jm CO; band generated via 
irradiation of CO; ice thin film with 100 keV H* protons, scaled in optical 
depth units (—In(R/Ro), where R is the reflectance of the film and gold 
substrate and Ro is the reflectance of the bare gold substrate; Raut & 
Baragiola 2013). All spectra are normalized to 1 at their maxima. The offset 
between the CO; features and Ariel’s 4.90 jm band supports the interpretation 
that this feature results from collective oscillations across a crystalline CO; ice 
lattice (Bini et al. 1991; Figure 2). It is possible that radiolytically produced 
CO; contributes to the broad short-wavelength side of Ariel’s 4.90 jum band. 
However, additionally experiments are required to determine whether this 
broadening of Ariel's 4.90 jum band results from mixing with CO; or other 
variables such as the temperature of CO; on Ariel. 


guest molecules in clathrates, however, would be more resistant 
to these processes and might persist. CO clathrates exhibit a 
"CO v3 mode near 4.685 um (Dartois 2011), potentially 
matching a sideband on Ariel's 4.67 um feature (Figure 3). 
Furthermore, the detection of CO in Enceladus's plume 
material (e.g., Peter et al. 2024) demonstrates that CO can 
survive in warm icy satellite interiors. 


5. Summary and Conclusions 


We analyzed JWST/NIRSpec spectral observations col- 
lected over Ariel's leading and trailing hemispheres, revealing 
the presence of crystalline CO; ice and CO ice mixed with 
CO». We detected a suite of other spectral features that might 
result from carbonates, C4O», and nitriles. The detected CO, 
ice is likely concentrated in thick deposits, possibly mixed with 
a small amount of CO, but the physical state of these deposits 
and whether they are associated with geologic conduits to 
Ariel’s interior are difficult to determine with these disk- 
integrated spectra. Spatially resolved NIR spectra collected by 
an orbiter making close flybys of the Uranian moons are 
required to explore the spatial relationship between volatiles 
and geologic features and confirm whether some carbon oxides 
originate in Ariel’s interior (e.g., Beddingfield et al. 2020; 
Cartwright et al. 2021; Leonard et al. 2021; Cohen et al. 2022; 
National Academies of Sciences, Engineering & 
Medicine 2023). 


Cartwright et al. 


We found no compelling evidence for NH-bearing species, 
hydrocarbons, or H5O5 in the G395M data, although NH- 
bearing species could be present but obscured by strong 3 wm 
H50 ice absorption. The lack of hydrocarbon features suggests 
that amorphous carbon dominates the low albedo material that 
is well mixed with H5O ice in the Uranian moons’ regoliths, 
possibly also mixed with some "amorphous silicates" (Cart- 
wright et al. 2018). The apparent absence of H5O» suggests that 
the irradiation environment at Ariel may be fairly quiescent, or 
extensive CO, deposits limit interactions between charged 
particles and underlying H20 ice. 
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Appendix 
A.1. Data and Methods: NIRSpec Observations 


Data processing utilized the Science Calibration Pipeline 
v1.13.4 with CRDS context jwst 1214.pmap to process raw 
uncal data into 53d spectral cubes for each of the four dithers 
(Bushouse et al. 2023). The pipeline was run using the default 
parameters and including the NSClean routine to remove the 1/ 
f pattern noise (Rauscher 2024). The spectral extraction used a 
"template PSF-fiting" routine. The wavelength grid was first 
computed using the CRVAL3 and CRDELT3 header keywords, 
with each wavelength corresponding to a specific slice in the 
data cube. An initial by-eye estimate was made for the centroid 
position of the target (Ariel is a point source given the 071 
NIRSpec pixels). The background was then calculated as the 
median of all pixels 75 pixels from the centroid and subtracted 
from all pixels in the slice. The "template PSF" was then 
calculated by taking the median of a moving 21-slice window. 
A 9x9 pixel box was cut out around the centroid of the 
template PSF, normalized to unity within the box, and 
iteratively fit to the slice in the middle of the 21-slice window 
using the scipy.optimize.minimize function and the Nelder- 
Mead algorithm. The two fit parameters were the flux scaling 
factor and the background, which were used to construct the 
best-fit model. The 1D spectrum was constructed by extracting 
the flux within a 3.5-pixel-radius circular aperture, centered on 
the centroid in each slice. To remove the solar component, the 
four dithers were medianed and then divided by a median 
G395M spectrum of P330E, a well-established spectrophoto- 
metric calibration star (GOV, Vmag 13.028 + 0.004; e.g., 
Bohlin & Landolt 2015). The individual P330E spectra were 
computed using the same template PSF-fitting routine 
described above. Uncertainties for the Ariel and P330E spectra 
were computed as the median absolute deviation within each 
wavelength bin, with the uncertainties propagated to the final 
P330E-divided spectra (Figure 1). 
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A.2. Data and Methods: Radiative Transfer Modeling 


The synthetic spectra reported in this study were generated 
by adopting Mie scattering theory (e.g. Bohren & Huff- 
man 1983) to calculate the single scattering albedo (Wo) for 
each component using their real, n, and imaginary, k, parts of 
the complex refractive index, derived from laboratory experi- 
ments for crystalline CO» ice (150K, Hansen 1997; 70K, 
Gerakines & Hudson 2020; 21K, Quirico & 
Schmitt 1997a, 1997b), CO ice (Robert Brown, private 
communication; 20 K, Gerakines et al. 2023), crystalline H2O 
ice (80 K; Mastrapa et al. 2009), and amorphous carbon (room 
temperature; Rouleau & Martin 1991). These c values were 
then passed to Hapke equations that calculate geometric albedo 
as a function of wavelength (Hapke 2012). The program 
generates one-layer models of end-member species (like those 
shown in Figure 2) or particulate or areal mixtures of various 
components (Figure A1). Minor resonances in synthetic spectra 
generated using Mie-derived wọ are addressed by calculating 
albedo using a range of grain sizes (410% spread in 
diameters), which are then averaged at each wavelength step 
in the final model. Although Mie scattering theory only 
approximates the structure of planetary regoliths, it is widely 
used in radiative transfer models to simulate the surfaces of icy 
bodies. More details and caveats on this program and its prior 
application to simulate the spectral properties of the Uranian 
moons' surfaces are provided in, e.g., Cartwright et al. (2023). 


A.3. Results and Analyses: Comparison between Ariel's 
4.90 um Band and Radiolytically Formed CO; 


Here we report laboratory spectra of radiolytically generated 
CO; compared to Ariel’s 4.90 jum band (Figure A2). 


A.4. Results and Analyses: Ariel's CO» Scattering Peak 
Compared to the Indices of Refraction for CO; 


Here we report an arbitrarily scaled NIRSpec spectrum of 
Ariel’s trailing hemisphere compared to the indices of 
refraction for CO , measured by Gerakines & Hudson (2020; 
Figure A3). 


A.5. Results and Analyses: Ariel Trailing/Ariel Leading 
Spectral Ratio 


Here we report a ratio between the spectra collected over 
Ariel's trailing and leading hemispheres (Figure A4). 


A.6. Results and Analyses: H505 and Ariel's 3.6 um Feature 


To investigate whether HO, is present on Ariel and to 
measure its 3.6 um feature, we fitted fourth-order polynomial 
models to the spectra of Ariel's leading and trailing hemi- 
spheres between 3.45 and 3.7 um and measured the resulting 
continuum-subtracted features (Figure A5). We find no 
evidence for the 3.505 um H202 feature detected on Europa 
(e.g., Carlson et al. 1999), Enceladus (Newman et al. 2007), 
and Charon (Protopapa et al. 2023, 2024; Figure A5). The 
continuum-subtracted data show weak features centered near 
3.55 and 3.60 um on Ariel’s trailing hemisphere. Whether these 
two features result from different species or are two subtle 
lobes of the same component is difficult to discern. We 
consider candidate species for the 3.55 and 3.60 jum features in 
Appendix A.7. 
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Figure A3. The real “n” (blue) and imaginary “k” (red) complex refractive 
indices for crystalline CO; ice measured in the laboratory at 70 K (Gerakines & 
Hudson 2020) compared to arbitrarily scaled spectra (1o errors) of Ariel's 
leading and trailing (offset upward for clarity) hemispheres. Dashed lines 
indicate spectral features identified on Ariel that result from CO, and "CO; 
ice and correspond to changes in n and k as a function of wavelength. Ariel's 
?CO, scattering peak starts close to 4.175 jum, where k substantially increases 
above n, and ends near 4.265 jum, where k drops off steeply and n increases 
above 6. The 4.30 jum 60? C 50 and 4.38 jum CO, features coincide with 
small increases in k above zero and slight dips in n. 


A.7. Discussion: Candidate Constituents 


Of the 15 features we identified and measured, seven are 
confidently attributed to H5O, '*CO>, or '*CO (Table 1). The 
other eight are centered near 3.33, 4.02, 4.15, 4.30, 4.41, 4.47, 
4.59, and 4.78 um. The 3.33 and 4.15 um features are likely 
dominated by CO», but the 3.33 um feature may include 
species exhibiting C-H stretching modes (e.g., Clark et al. 
2009; Grundy et al. 2002; Dartois et al. 2010), whereas the 
4.15 um feature may also include contributions from deuterated 
water ice (D;O/HDO); Clark et al. 2019). 

Ariel’s 4.02 um feature might result from a broad v, + v4 
combination mode expressed by CO3-bearing minerals (Hex- 
ter 1958; Bishop et al. 2021). Additionally, S-bearing species 
have been suggested to explain Callisto's broad 4.02 um 
feature, including disulfanide (HS5) and hydrogen disulfide 
(H2S2; Cartwright et al. 2020c), or perhaps even a wavelength- 
shifted sulfur dioxide (SO2) band (e.g., McCord et al. 1998). 
H,CO3, formed via irradiation of H5O mixed with COs, 
expresses a broad band between 3.8 and 3.9 jum due to the 
CO; anion (e.g., Hage et al. 1998) that might contribute to 
Ariel's 4.02 jum, assuming that some process is able to shift its 
band center closer to 4 jum. Conversely, Ariel’s 4.02 um feature 
might be associated with the Christiansen band exhibited by 
"CO: ice near 4.07 um (seen in the CO» spectral models 1—3 
shown in Figures 2 and Al). However, there is a sizable 
wavelength gap (70.05 jum) between the Christiansen features 
displayed by these spectral models and Ariel's 4.02 um band, 
and it is uncertain why the Christiansen band would exhibit 
such a large shift, especially given that the feature is not shifted 
in spectra of CO»-rich TNOs (Figure A6). 
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Figure A4. Ratio between Ariel’s trailing and leading hemispheres, normalized to 1 at 4 jum. Dashed lines highlight the spectral features identified in the Ariel data 
(Figure 1) that show notable trailing/leading hemispherical asymmetries, including subtle features near 3.01, 3.33, 4.84, and 4.93 jum that could result from deposits of 
CO; ice upward of 10 mm thick (Section 4.1, Figure 4). Candidate species are described in Appendix A.7. 


Because of the strength of H5O ice absorption and the shape 
of the H,O ice continuum in the 4 um wavelength range, exact 
identification of Ariel's 4.02 um band center is difficult, and it 
may be centered at slightly shorter wavelengths, between 3.98 
and 4 um. Similarly, Ceres exhibits a 4 jum feature that shifts 
between 3.95 and 4.02 um, with Mg-Ca carbonates shifting its 
4m feature to the global average of 3.95 m and Na 
carbonates shifting the band feature to 4.02 um (e.g., Carrozzo 
et al. 2018). Therefore, if Ariel's 4.02 jum band is in reality 
centered at slightly shorter wavelengths, it would still be 
consistent with the broad v, + v4 combination mode expressed 
by carbonate minerals. 

Ariel's 4.30 um feature is embedded on the long-wavelength 
end of its 4.27 um "CO, ice band and probably results from 
the CO, isotopologue '°O'*C'8O0 (Bennett et al. 2010b). 
Although amorphous CO; could hypothetically contribute to 
Ariel's 4.30 jum band, it is not expected to be present at Ariel's 
peak surface temperatures (e.g., Escribano et al. 2013). 

The 4.41 and 4.47 um features could result from irradiation 
of CO and CO, mixtures, forming C30, and other higher-order 
carbon chain oxides (e.g., Gerakines & Moore 2001; Strazzulla 
et al. 2007). Alternatively, these two features could result from 
irradiation of carbon oxides and ammonia (NH3), forming CN- 
bearing compounds (e.g., Strazzulla et al. 2007). The 4.41 uum 
feature might also (in part) result from 16936159 ice (Bennett 
et al. 2010b). 

Ariel's 4.59 um feature could result from irradiation of H20, 
CO», and NH}, forming OCN , a key tracer of nitrogen- 
bearing species in the interstellar medium and protoplanetary 
disks (e.g., McClure et al. 2023). If OCN is confirmed, then a 
complementary NH4 feature should be present near 6.85 jum 
(e.g., Grim et al. 1989; Palumbo et al. 2000; Bennett et al. 
20102). 

Ariel's 4.78 um band most likely results from 13CO (e.g.. 
Bennett et al. 2010b), given that it is only detected on Ariel's 
trailing side, where the '2CO feature is dramatically stronger 
(see Appendix A.8 for discussion of CD / ?CcO isotopic 
ratios). Nonetheless, the v3 mode of hydrogen cyanide ice 
(HCN; e.g., Gerakines et al. 2022) might contribute to this 
feature as well, if nitriles are present. 

Subtle features near 3.55 and 3.60 um on Ariel’s trailing side 
(Figure A5) may result from irradiation of CO; and H50, 
possibly generating formaldehyde (H2CO; M. Moore et al. 
2004), formic acid (CH5O»; Bisschop et al. 2007), or perhaps 
oxalates (C204; Applin et al. 2016). A 3.6 um feature on 
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Saturn’s moon Phoebe has also been attributed to HDO/D5O 
ice (Clark et al. 2019), and perhaps deuterated water ice 
contributes to Ariel's 3.6 jum feature as well. 

Finally, we identified, but did not measure, three other very 
subtle features centered near 3.01, 4.84, and 4.93 jum on Ariel's 
trailing hemisphere, which probably result from CO, (Figures 4 
and A4; Hansen 1997; Quirico & Schmitt 19972). 

Disentangling these various components requires follow-up 
laboratory experiments and associated modeling efforts. In 
particular, there is a need to determine whether NH, on Ariel 
may have been radiolytically converted into nitriles, thereby 
explaining the absence of N-H stretching modes in the 
NIRSpec data reported here. Additionally, Ariel’s prominent 
double-lobed scattering peaks near 4.20 and 4.25 jum share 
morphological similarities to double-lobed CO» gas emission 
peaks detected at Callisto (Carlson 1999) and, more recently, 
the Centaur 39P/Oterma (Pinto et al. 2023). However, there is 
a notable wavelength shift between these CO» gas emission 
peaks (~4.2-4.3 um) and Ariel’s scattering peaks 
(4.16-4.26 jum), which is more consistent with scattering 
within a thick layer of CO; ice (Model 3 in Figure 2; 
Figure A3). 


A.6. Discussion: Comparison between Ariel and CO;-bearing 
TNOs 


Here we compare NIRSpec data of Ariel's leading and 
trailing hemispheres to spectra of three TNOs observed as part 
of GO program 2418 (Pinilla-Alonso et al. 2021; Figure A6). 
These TNOs are representative of the three spectral types 
identified by this program (de Pra et al. 2024). 


A.9. Discussion: '^CO/"CO Isotopic Ratios 


Prior studies have used remotely sensed H2O, HDO/D5;O 
ice, PCH, 13 CH4, 12CH;D, i CO», and ?CO, spectral features 
to estimate D/H and '^C/P^C isotopic ratios and gain insight 
into the formation conditions for different icy bodies and 
possible endogenic sources of material on their surfaces (Clark 
et al. 2019; Cartwright et al. 2024; Glein et al. 2024; Grundy 
et al. 2024b). Ariel's strong CO, scattering peaks make 
reliable determination of its 8CO, / “CO; isotopic ratios 
difficult from analysis of v; band parameters for '“CO3, and 
we did not attempt to retrieve a CO, abundance from other 
spectral features (e.g., at 4.90 um). However, using the 
4.67 um '*CO and 4.78 um (likely) ^CO bands detected on 
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Figure A5. (a) Spectra of Ariel’s leading (purple) and trailing (red) hemispheres, scaled to arbitrary reflectance, fit by fourth-order polynomials simulating their 
continua. Dashed lines highlight the canonical central wavelength position of H2O, measured in the laboratory (3.505 um) and two other subtle features centered near 
3.545 and 3.602 jum we have identified on Ariel's trailing hemisphere. (b) Continuum-subtracted spectra of Ariel’s leading (purple) and trailing (red, offset upward by 
0.0025) hemispheres, spanning the 3.4—3.7 jum wavelength range. A two-Gaussian model (black) has been fit to the weak 3.55 and 3.60 jum features observed on 
Ariel's trailing hemisphere that are not seen on its leading hemisphere. Neither spectrum exhibits a 3.505 jum band resulting from a H2O, combination mode (Bain & 


Giguére 1955). 


Ariel's trailing hemisphere, we calculated "cuo "CO band 
area and depth ratios of 0.045 + 0.010 and 0.146 + 0.023, 
respectively. These results are comparable to Phoebe's 
100,/?CO, band area (0.059+0.005) and depth 
(0.186 + 0.036) ratios (Cartwright et al. 2024), and the ratios 
on Ariel’s leading hemisphere appear qualitatively to be higher, 
hinting that Ariel may be significantly enhanced in "°C relative 
to "terrestrial" values (Pc/ 12C ~ 0.010-0.012) exhibited by 
most bodies in the inner solar system and Saturn's rings and 
regular satellites (e.g., Clark et al. 2019). Whether "C/'^C 
ratios derived from CO at Ariel would yield similar results to 
carbon isotopic ratios derived from CO, at Phoebe is uncertain. 
Furthermore, CO is likely contaminated by CO, at Ariel, and 
its 4.78 wm band may be enhanced by HCN ice, assuming that 
nitriles are present, complicating analysis of its carbon isotopes. 
We report a preliminary analysis of possible "°C enrichment of 
CO and CO, ices on Ariel in Appendix A.10. More detailed 
modeling of Ariel’s carbon isotopes and their spectral 
expression at infrared wavelengths is likely required to 
corroborate the CO band ratios present here, to convert band 
ratios to number ratios, and to fully explore their significance, 
e.g., in terms of formation of its carbon source material (e.g., 
Cartwright et al. 2024). 


A.10. Discussion: Modeling CO, Enrichment 


The pronounced 4.38 um absorption due to CO; and high 
apparent ratios of PCO / CO (assuming that CO is a radiolytic 
product of CO2) provide hints that Ariel’s surface CO» ice is 
enriched in "C, perhaps substantially. Here we seek to outline 
explanations of how CO» could acquire a large C enrichment 
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via a balance between production and escape processes. We 
then consider a different scenario in which there is long-term 
isotopic fractionation of CO; driven only by escape. 

Since CO) is concentrated on the trailing hemisphere, it may 
be produced by radiolytically driven oxidation reactions of, 
e.g., organic compounds. An alternative or complementary 
source of CO, would be outgassing from the interior (e.g., a 
soda ocean). It may not be essential to specify the nature of the 
CO, source since both options would provide CO, with 
roughly similar C isotope ratios (see below). The exception 
would be if radiolytic processes result in isotopic fractionation, 
but high-energy processes are not generally associated with 
strong isotope effects because they usually lack selectivity. We 
can formulate two equations for the evolution of CO; 
isotopologues in the surface environment of Ariel: 


dNa4 


= Ja — ka [PCOo]exob (Al) 
dt 
and 
dN. 
= Jis — Kas[ PCO» exon, (A2) 


where N represents the number of molecules, subscripts refer to 
the masses of CO; isotopologues, t stands for time, J designates 
production rate, k represents a rate constant for atmospheric 
escape, and brackets indicate the number density of the CO; 
isotopologue of interest at the exobase. We adopt the classic 
concept of the exobase as a dividing line between nominally 
collisional and collisionless parts of the atmosphere. We wish 
to explore the simplest case in which a steady state could be 
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Figure A6. Comparison between NIRSpec IFU spectra and 1c uncertainties for Ariel’s leading (purple) and trailing (red) hemispheres and the TNOs 47171 Lempo 
(1999 TC36) (1), 2013 LU28 (2), and 2004 XA192 (3), selected to represent members of the “bowl,” “cliff,” and “double-dip” TNO types (Pinilla-Alonso et al. 2021; 
de Pra et al. 2024), respectively. All spectra are normalized to 1 at 4.15 jum and offset vertically for clarity. Similar to other double-dip TNOs, 2004 XA192 exhibits 
strong CO; and CO, bands, with a CO, scattering peak similar to Ariel's 4.202 um peak. However, none of the TNOs display a strong 4.252 jum peak, unlike 
Ariel. 2004 XA192 also exhibits a strong ~4.07 um band resulting from the Christiansen effect, where n = 1 and k = 0 in crystalline CO» ice (Figure A3). This 
Christiansen band is notably absent from the Ariel spectra and synthetic spectra composed of large grain sizes (diameters ~100 jum; Gerakines & Hudson 2020), 


generated using our Hapke—Mie approach (Appendix A.2). 
established. Equations (A1) and (A2) then imply that 


"CO; _ Jas/Ja _ Rpro 
12CO Jop — Kas/Kaà Qese 


(A3) 


where Ry, stands for a ratio characteristic of a given 
production process and o&,. is a fractionation factor due to 
CO; loss by atmospheric escape. 

We treat Ry, as a fixed parameter that may assume =e / 12c 
ratios between those measured in primitive organic matter (1 / 
92; Alexander et al. 2007) and carbonate minerals (1/83; 
Fujiya et al. 2023) in carbonaceous chondrites. The former 
would reflect an organic source of Ariel's CO; ice (e.g., 
carbonaceous dust), while the latter can serve as an analog for 
carbonate-bearing subsurface fluids that might erupt CO». The 
second value would be inherited from protosolar nebula CO2 
ice that may have been accreted by Ariel at the time of its 
formation. The overall range of C / 12C ratios considered here 
encompasses a broad range of primordial carbon isotope ratios 
observed across the solar system. 

We consider CO; loss by Jeans escape. This mechanism has 
the potential to discriminate between species with different 
masses (hydrodynamic escape and sputtering are other 
possibilities, but they are less effective than Jeans escape at 
fractionating isotopes). For Jeans escape, 


44 É + x 


Qese = 


ex (X = A’). (A4) 
alt 4 ur P (444 45 
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where A®%® represents the escape parameter of species i at the 


exobase, as defined below: 
GMmij 


exob 
A = , 
k B Texob Texob 


(A5) 


where GM = 8.23 x 109 mê s? m; corresponds to the mass of 
a molecule of species i; kg = 1.381 x 10? kgm?s ?K 5; 
Texop denotes the exobase temperature, which is assumed to be 
the same as the rest of the atmosphere (as a first-order 
approximation); and Fexob indicates the distance between the 
exobase and the center of Ariel. 

The isotope ratio at the surface can be related to that at the 
exobase via 


BCO,) | 1 (PCO 
DOO, Jut — Odit \7CO2 Jaron | 


where o; denotes an isotope fractionation factor due to 
molecular diffusion between the surface and exobase. We 
assume that turbulent mixing is unimportant in Ariel's tenuous 
atmosphere (i.e., the homopause is always at the surface). For 
an isothermal atmosphere, the diffusive fractionation factor can 
be computed using 


(A6) 


GM (mas — maa) raus — rout) 
kg Tim 


Qai = ew } (A7) 


where F uf =579 x 10°m. We introduce Tam aS a mean 
temperature at altitudes between the surface and exobase. 
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We do not know the location of Ariel’s exobase. One 
possibility is that Ariel has a surface-bound exosphere. We can 
determine whether this is likely by computing the surface 
pressure (Pourr) that would satisfy the following relation: 


2 2 
deo, sut Burt 


=> 1, 
GMmco, 


(A8) 
where dco, corresponds to the kinetic diameter of CO; 
(4.53 x 10? m; Haynes 2016). The pressure would need to be 
«2.8 x 1078 Pa for the exobase to be at the surface. If there is a 
sufficient CO, supply (as suggested by the prominent CO, 
features in the JWST data), then we may assume vapor pressure 
equilibrium to set a constraint on the maximum temperature 
that permits a surface-bound exosphere. We use the equation 
from Fray & Schmitt (2009): 


5 


In pco, sa (bar) = Ao + » Aj Turt (KY, 
j=l 


(A9) 


where pco»,4 designates the saturation pressure of CO, ice, 
Tsuf indicates the surface temperature, Ao = 14.76, 
A, =—2571K, A55 —7.781 x 10* K^, A4—4.325 x 10°K’, 
A4 — —1.207 x 108 K^, and A5— 1.35 x 10? K". These para- 
meters are applicable to temperatures between 40 and 195 K 
(Fray & Schmitt 2009). By calculating the saturation pressure 
as a function of temperature, we find that the left-hand side of 
Equation (A8) is «1 (meaning that the mean free path is longer 
than the scale height) when T; « 71 K. Hence, the exobase 
will be at the surface (i.e., the whole atmosphere is an 
exosphere) when the temperature is below 71 K; in this case, 
Qaire = 1, and we can solve Equations (A3)-(A6) to determine 
the O0, / "co, ratio at the surface of Ariel. 

At higher surface temperatures, the exobase will be at higher 
altitudes. It is not straightforward to simulate the structure of 
this type of atmosphere; this is an active area of research 
(Mogan et al. 2020). However, detailed modeling is probably 
not needed for our application since we lack specific values of 
the 13CO, / 12C0, ratio on Ariel to explain (see Appendix A.8). 
Instead, we have the more basic goal of exploring the 
sensitivity of the isotope ratio to unknowns that influence 
isotopic fractionation. The two parameters of interest here are 
the exobase altitude and temperature of the atmosphere. In the 
present treatment, we still assume that the atmosphere is 
isothermal, but it could have a different temperature from that 
of the surface. Parameters rexop and Tatim in Equation (A7) can 
be varied to assess how much they could change ovi. 
Equations (A3)-(A6) can then be used to estimate the steady- 
state P? C/ 12C ratio of Ariel’s surface CO, ice at Touf > 71 K. 
We do not account for any carbon isotopic fractionation 
between CO2 gas and ice at the surface, as laboratory 
experiments show that there is no more than ~0.4%o 
fractionation down to 130 K (Eiler et al. 2000). 

As shown in Figure A7, CO; derived from organics or 
carbonates can be significantly enriched in ^C in the surface 
environment of Ariel due to atmospheric escape. The predicted 
enrichment is about 30%, i.e., 300%c. This is very large. For 
comparison, biological production of CH, from CO, is 
generally considered to have a large isotope effect, with a 
fractionation that can reach ~80%o (Milkov & Etiope 2018). 
Our predictions appear to be qualitatively consistent with the 
ease of finding features from CO, and PCO in the JWST 
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Figure A7. Carbon isotopic composition (light/heavy) of steady-state CO» at 
the surface of Ariel. These results show how input values of the ^C/'^C ratio 
from organics or carbonates (dark horizontal lines) are shifted downward by 
Jeans escape. For each set of lighter-colored curves, we consider a case with a 
surface-bound exosphere at all temperatures and Tam = Tsurf (dotted), a case 
where the exobase altitude has been increased to 2000 km while keeping 
Tam = Tsur¢ (dashed), and a case where the exobase altitude remains at 2000 
km but the atmospheric temperature is 20 K lower than the surface temperature. 
The offsets for these cases are intended to be illustrative, and other values are 
possible. 


data. It can be seen in Figure A7 that increasing the exobase 
altitude may not have a significant effect on the isotope ratio, at 
least up to 2000 km. The atmospheric temperature can have a 
larger effect in enriching CO, at the surface. Yet its effect on 
diffusive fractionation seems less important than Jeans 
fractionation, which is responsible for most of the total 
fractionation in our test cases. Isotopic fractionation by 
diffusion is limited by large scale heights of Ariel’s 
atmosphere. 

The spurious drops at 71 K in Figure A7 arise from the 
arbitrary assumption that the exobase is at the surface at 
temperatures below 71 K and at an altitude of 2000 km at 
higher temperatures. Although assuming a constant exobase 
level isolates its effects on the '^C/'?C ratio, in a more realistic 
model the exobase would progressively increase in altitude 
from the surface level with increasing surface temperature and 
pressure. 

Whether a steady state is reached is an open question. In 
another end-member scenario, one can assume progressive 
depletion of a primordial surface CO, inventory, as a fraction 
sublimates each summer to form gaseous CO» that undergoes 
isotopic fractionation due to Jeans escape. From Equation (A4), 
Qesc % 0.72 assuming Taf —55K and a surface-bound 
exobase, or 0.78 at Tsuf = 71 K. Assuming escape of a modest 
fraction fes: — 0.1 ppm of all the surface CO, each Uranian 
year, the ^C / 12C ratio of primordial CO; would increase by a 
factor of (fesc/ Qese) (4.57 x 10? /86) = 1 (i.e., 7000960) over the 
4.57 Gyr of solar system history, assuming a constant orbital 
period for Uranus of 86 Earth years. This requires CO; to be 
gaseous for part, and only part, of the Uranian year. 

This calculated enrichment in ^C is even more considerable 
than in the steady-state scenario, as the ever-increasing "C / PO 
ratio of the surface CO; is not diluted by constantly 
replenishing CO, with a canonical ?C/'^Cz: 1/90. This 
enrichment is of the same order of magnitude as the factors 
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of ~4.5 and 15 inferred, respectively, from the band area and 
depth ratios of '^CO and '*CO on Ariel (see Appendix A.9). 
The corresponding Fese o ppm is consistent with a CO, 
initial inventory (1 + fa 6 ~ 200 times greater than the 
present-day surface abundance; i.e., equivalent layer thick- 
nesses on the order of 10m or larger based on the >10 mm 
present-day deposit thicknesses discussed in Section 4.1. 
Compared to water inventories equivalent to global ice shell 
thicknesses of 100km or more (e.g., Castillo-Rogez et al. 
2023), this appears compatible with bulk CO5/H5O abundance 
ratios of less than a few mol% based on observations of comets 
(Mumma & Charnley 2011), which may be indicative of the 
composition of icy material accreted by Ariel. Such bulk 
compositions allow for greater values of fesc, which would be 
required to obtain the same "^C enrichment in intermediate 
scenarios involving both primordial inventory depletion and 
some degree of endogenic (and, for CO, radiolytic) replenish- 
ment. We leave to future work the questions of what the 
absolute rate of CO; escape might have been through time, 
whether it is consistent with the removal of almost all of Ariel’s 
surface CO, inventory, and how much radiolytic production 
may be implied to account for trailing versus leading 
hemisphere differences in the apparent abundance of CO;. 
Future constraints on the !50/!^O ratio may enable further 
discrimination between steady-state and limited replenishment 
scenarios, as the larger mass difference between '°OC'*O and 
!'60C!'60 can be expected to magnify the degree of 
fractionation. 

The models outlined above may have more general 
applicability to other outer solar system bodies with abundant 
‘SCO . In addition to the phenomena considered in this 
appendix, it may be useful to consider how lateral hetero- 
geneities in the availability of CO» ice, seasonal effects on 
surface temperatures, kinetic inhibition of CO» sublimation, 
and possible vapor pressure isotope effects at low temperatures 
might come together to create a rich isotope geochemistry that 
we are just beginning to glimpse. 
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